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Abstract 

The  mechanical  energy  balance  for  the  sea  ice  cover  of  the 
Beaufort  Sea  portion  of  the  Arctic  Ocean  has  been  used  to  evaluate  a 
parameter  study  on  the  strength  of  a  plastic  sea  ice  model.  The  results 
of  this  study  show  that  a  hardening  plastic  model  with  a  diamond  shaped 
yield  surface  and  a  compressive  strength  of  1.0  to  1.5  x  10^  Nm  1  will 
simulate  the  winter  behavior  of  ice  in  the  Beaufort  Sea  accurately. 

The  mechanical  energy  (power)  budget  utilizing  a  hardening  plastic 
model  was  examined  for  a  seventeen  day  period  during  the  winter  of  1976. 
The  study  period  included  two  storms  (extreme  ice  motion  events)  during 
which  the  power  from  the  atmosphere  is  transferred  to  the  sea  ice  cover 
and  redistributed  horizontally  by  the  divergence  of  the  stress  flux, 
concentrating  the  dissipation  of  mechanical  energy  in  the  sea  ice 
cover  along  the  north  slope  of  Alaska.  These  results  show  how  the  power 
from  the  atmosphere  over  the  entire  Beaufort  Sea  can  be  concentrated  to 
ridge  and  raft  sea  ice  and  generate  acoustic  energy  in  the  ocean  along 
the  margins  of  the  Beaufort  Sea. 
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Introduction 

In  temperate  regions  the  transfer  of  mechanical  energy  between  the 
atmosphere  and  the  ocean  has  been  studied  (e.g.,  Holland,  1978).  At 
high  latitudes  where  the  ocean  is  ice  covered  during  some  seasons  of  the 
year,  the  transfer  of  mechanical  energy  is  altered  and  only  recently 
(Coon  and  Pritchard,  1979)  has  the  concept  of  mechanical  energy  balance 
been  introduced  in  such  a  way  as  to  include  sea  ice.  This  concept 
provides  a  simple  description  of  the  transfer  of  energy  through  the 
atmosphere/ice/ocean  system  and  explains  how  the  ice  cover  dissipates 
energy  by  deformation  and  transfers  energy  horizontally  by  stress  flux 
divergence. 

The  dissipation  of  mechanical  energy  in  the  sea  ice  cover  results 
in  ridging,  rafting  and  the  generation  of  acoustic  energy.  Coastal 
processes,  such  as  beach  erosion  and  sediment  transport,  are  also 
altered  by  the  presence  of  the  ice  cover  and  its  effect  on  the  mechani¬ 
cal  energy  budget.  Therefore,  understanding  how  the  mechanical  energy 
from  the  atmosphere  is  transferred  to  the  sea  ice  cover  and  subsequently 
redistributed  and  dissipated  within  the  ice  cover  leads  to  an  under¬ 
standing  of  these  physical  phenomena. 

In  this  paper  a  seventeen  day  period  during  the  winter  of  1976  is 
examined  for  the  Beaufort  Sea.  The  energy  budget  is  used  to  evaluate  a 
parameter  study  of  strength  and  yield  surface  shapes  for  a  plastic  ice 
sea  model.  Material  parameters  in  the  model  are  optimized  and  the 
tuned  model  is  subsequently  used  to  study  the  energy  budget  of  the 
atmosphere/ice/ocean  system.  This  study  shows  how  the  energy  from  the 
atmosphere  distributed  over  the  ice  cover  of  the  Beaufort  Sea  can  be 
concentrated  and  dissipated  along  the  north  slope  of  Alaska  and  the 
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Canadian  Archipelago.  During  these  periods  of  high  energy  dissipation, 
severe  ridging,  rafting  occur  in  these  marginal  seas. 

Test  Conditions 

From  January  27  through  February  12,  1976,  winds  in  the  Beaufort 
Sea  (see  Figure  1  for  the  region  considered)  cause  ice  behavior  that  is 
both  interesting  and  useful  for  testing  ice  model  response.  During  this 
time  period  two  storms  pass  through  the  Beaufort  Sea  causing  large 
motions,  first  to  the  west  and  later  to  the  east.  The  second  storm 
lasts  beyond  the  test  period.  During  that  time  the  entire  ice  cover 
remains  motionless,  even  though  strong  winds  persist.  A  brief  des¬ 
cription  of  the  day-to-day  activity  follows. 

On  January  27  and  28,  while  moderate  easterly  winds  blow  in  a  band 
about  200  km  wide  off  the  Alaska  north  slope,  the  ice  remains  stationary. 
On  January  29,  a  high  pressure  region  begins  to  move  eastward  across  the 
northern  edge  of  the  Beaufort  Sea  (above  the  grid  boundary  Figure  1). 
Westward  ice  motions  begin  at  the  western  edge  of  the  region.  By 
January  31,  the  entire  region  is  moving  westward  (from  Barrow  to  Banks 
Island)  except  for  a  wide  belt  of  fast  ice  along  the  Alaska  north  slope. 
The  fast  ice  region  is  shoreward  of  the  flaw  lead  running  from  Barrow  to 
Banks  Island  and  about  130  km  offshore  near  Prudhoe  Bay.  During  this 
time  a  series  of  north-south  running  leads  opens  uniaxially  with  the 
lead  system  appearing  first  at  the  western  edge  of  the  region  and  then 
moving  eastward  with  the  storm.  Finally  a  50  km  wide  lead  appears 
along  Banks  Island.  The  winds  die  down  on  February  1,  and  come  about. 
There  is  little  activity  through  February  5. 

On  February  6,  a  low  pressure  area  to  the  north  of  the  region  of 
interest  causes  westerly  winds  to  rise  over  most  of  the  region.  The  ice 
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begins  to  move  eastward.  Then  winds  strengthen  to  give  an  air  stress  of 
about  0.4  Pa  and  the  ice  moves  toward  Banks  Island  with  substantial 
shearing  along  the  Alaska  north  slope.  The  previous  fast  ice  region 
then  moves  with  the  pack.  By  February  8,  the  ice  nearest  Banks  Island 
stops  moving  because  the  open  water  and  thin  ice  have  been  eliminated 
by  ridging.  Ice  strength  increases  so  that  the  stress  can  resist  the 
wind  force.  This  region  spreads  so  that  by  February  10,  the  entire 
Beaufort  Sea  is  motionless  even  though  strong  winds  continued  to  blow 
after  the  end  of  the  test  period  on  February  12. 

During  this  time  period  there  are  twenty  data  buoys  and  three 
manned  camps  located  in  the  region  of  interest.  These  provide  accurate 
motions  and  winds  to  test  the  response  of  the  model  (Pritchard,  et  al., 
1977).  Therefore,  in  addition  to  the  energy  variables,  there  are  velocity, 
deformation  and  data  available  to  support  the  results. 

Sea  Ice  Model 

The  mathematical  model  under  consideration  in  this  work  is  a  modified 
version  of  the  AIDJEX  plastic  ice  model  (Coon,  et  al.,  1974;  Pritchard, 
1975;  Rothrock,  1975;  Thorndike,  1975).  In  the  momentum  balance  equa¬ 
tion  this  model  accounts  for  air  stress  T  ,  water  stress  T  ,  Coriolis 

~a  -w 

force  -mfkxv  ,  sea  surface  tilt  mgVH  and  internal  stress  divergence 
V*o  . 

mv  *  T  +  T  -  mfk*v  -  mgVH  (1) 

where  v  is  ice  velocity  and  m  is  area  mass  density.  The  ice  stress 
is  related  to  deformation  by  an  elastic-plastic  constitutive  law.  Yield 
strength  depends  on  the  instantaneous  thickness  distribution  at  each 
location  so  that  hardening  or  softening  can  occur  with  deformation.  The 
thickness  distribution  describes  the  relative  fraction  of  ice  of  each 
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thickness  category  and  the  concentration  of  each  category  changes  as  the 
ice  deforms  through  ridge  and  lead  formation.  Details  of  model  formula¬ 
tion  and  parameters  are  available  in  the  AIDJEX  final  report  to  the 
National  Science  Foundation  (Coon,  et  al.,  1978).  A  brief  description 
of  material  parameters  is  given  by  Coon  and  Pritchard  (1979).  In  the 
present  work  the  strength  parameter  variations  are  performed  for  a 
simpler  model  in  which  yield  strength  (magnitude  and  shape  of  yield 
surface)  is  assumed  fixed  for  a  given  calculation.  The  yield  surface  is 
then  varied  for  different  calculations  to  determine  the  effect  of  this 
parameter.  Previous  work  by  Pritchard  et  al.,  (1977)  has  shown  that  the 
ideal  plastic  model  (no  hardening/softening)  provides  daily  average 
velocity  and  strain  results  that  are  nearly  equal  to  those  obtained  for 
the  hardening/softening  model  if  a  sequence  of  daily  winds  is  used  in 
place  of  the  constantly  varying  winds.  The  shape  of  the  yield  surface 
will  be  varied  in  this  parameter  study  between  a  teardrop,  triangle  and 
diamond  shape  as  shown  in  Figure  2,  Boundary  layer  drag  coefficients 
are  unchanged  for  all  runs. 

Mechanical  Energy  Balance 

Consider  a  large  scale  element  of  ice,  say  tens  of  kilometers 

across,  composed  of  many  individual  ice  floes,  leads  and  ridges  (Figure 

3).  At  the  top  surface  of  the  element  of  ice,  the  wind  works  on  the  ice 
T 

at  the  rate  p  *  v  t  which  is  the  product  of  the  air  stress  acting 

a  ~  -a 

on  the  ice  and  the  ice  velocity.  Each  term  in  p^  may  be  considered 

an  average  over  the  area  which  is  reduced  to  a  point  by  a  standard 

limiting  argument.  At  the  bottom  surface  of  the  ice  element,  the  ice 

T 

works  on  the  ocean  at  the  rate  of  p  ■  -  v  T  where  T  is  the  traction 

o  -  ~w  ~w 

exerted  on  the  ice  by  the  ocean. 
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Around  the  sides  of  the  ice  element,  the  stress  field  produces  a 
traction  on  where  a  is  actually  a  stress  resultant  (the  Cauchy 
stress  in  excess  of  hydrostatic  equilibrium  integrated  through  the  ice 
thickness),  n  is  the  unit  outward  normal  to  the  side  of  the  ice  element. 

-s# 

The  rate  at  which  work  per  unit  area  is  done  on  the  ice  element  by  this 
traction  is 

Pf  =  J  f  yTon  d£.  (2) 

L 

where  L  is  the  curve  defining  the  horizontal  area  A  of  ice  element. 

From  the  Green-Gauss  theorem  this  power  input  at  the  sides  may  be 

T 

expressed  as  the  divergences  of  the  stress  flux  pf  *  V* (v  a)  ,  which 
is  the  horizontal  transfer  of  energy  by  the  nonhomogeneous  velocity  and 
stress  fields.  In  this  development,  the  velocity  field  is  assumed 
continuous. 

The  remaining  force  acting  on  the  ice  element  that  contributes  to 
the  mechanical  energy  budget  is  gravitational  and  is  important  only  when 
the  sea  surface  is  tilted.  The  rate  of  increase  of  gravitational  poten¬ 
tial  energy  as  the  ice  moves  up  the  sloped  surface  H  is 

T 

Pg  =  v  mgVH 

where  m  is  the  mass  of  the  ice  element  and  g  is  the  gravitational 
constant. 

The  mechanical  energy  balance  is  obtained  by  equating  the  resultant 

of  work  done  by  all  forces  acting  on  the  ice  element  to  the  sum  of  the 

rate  of  work  done  by  internal  stress  accompanied  by  deformation 

p^  *  tr  aD  where  D  is  the  stretching,  and  rate  of  change  of  kinetic 
1  T 

energy  K«^mvv  .  Therefore, 
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*  +  pd  -  ».  -  p„  -  pg  +  pf  <3> 

This  result  is  identical  to  that  of  Coon  and  Pritchard  (1979)  who  intro¬ 
duced  pw  *  pQ  +  Pg  as  the  total  power  transfered  into  the  ocean. 

Although  the  present  decomposition  into  power  input  to  the  ocean  by  the 
water  drag  force  and  the  rate  of  increase  of  gravitational  poten¬ 

tial  energy  is  useful  for  understanding  the  physical  importance  of  each 
term,  we  use  the  total  power  to  the  ocean  in  the  remainder  of  the 

analysis  primarily  because  p  is  small  and  there  is  no  difference  in 

S 

interpretation  of  results. 

In  the  present  work,  a  seventeen  day  period  is  studied  as  a  sequence 
of  daily  response.  For  time  resolution  on  the  order  of  one  day,  a 
quasi-steady  analysis  is  adequate.  Thus,  rate  of  change  of  kinetic 
energy  is  ignored  and  the  sequence  of  instantaneous  values  satisfies 

pa  +  pf  =  pw  +  pd  (4) 

The  power  input  to  the  ice  from  the  atmosphere  p  is  calculated 

a 

from  the  air  stress  T  =  p  c  I  I U  I  I B  U  where  U  is  the  geostrophic 

-a  a  a 1 1 ~g 1 1 ~a^g  -g  &  r 

wind,  B  is  a  rotation  operator  through  angle  ot(ccw)  and  pc  is 
air  density  times  a  drag  coefficient.  It  is  satisfactory  to  think  of 
the  air  stress  Ta  as  the  input  variable  (Pritchard,  et  al.,  1977). 

The  water  stress  is  T  =  p  c  I |v  -  v  I |B(v  -  v  )  where 

-W  W  WM~  ~g  '  1  ^  ^  -.g7 

3 

pc  *  5.5  kg/m  is  water  density  times  a  drag  coefficient  and  v  is  the 
w  w  ~  g 

geostrophic  ocean  current  that  is  assumed  to  be  constant  in  time  for  all 
calculations  (e.g.,  Pritchard,  et  al. ,  1977). 

The  global  form  of  mechanical  energy  balance  is 


P  +  Pr 
a  f 


P 

v 


+  pd 


(5) 


where  for  each  term 


P 


1 

A 


J  P  da 
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and  each  term  in  Equation  (4)  is  the  area  average  of  its  counterpart  in 
Equation  (3). 

One  additional  relationship  that  is  required  is  the  free  drift 
power  input  Pa£^  *  This  term  is  the  power  that  would  be  input  locally 
to  the  ice  from  the  atmosphere  if  the  ice  cover  were  present  but  ice 
stress  were  zero.  For  this  case  a  =  0  implies  p^  =  ^  0  so  that 

all  free  drift  energy  is  input  locally  to  the  ocean.  Thus, 


Pafd  Pwfd 


(6) 


and  a  similar  global  form  follows  directly. 


Results 

Data  from  ten  calculations  with  the  plastic  sea  ice  model  are 
presented  a)  to  show  how  the  energy  transfer  is  altered  by  changing 
strength  of  the  ice,  b)  to  show  how  the  information  may  be  used  to 
verify  and  tune  material  parameters,  and  c)  to  identify  important  con¬ 
tributions  to  the  mechanical  energy  budget  during  winter  when  ice  is 
compact  and  stress  is  important. 

Strength  Variations 

A  set  of  eight  simulations  during  the  first  eight  days  of  the  test 

period  (January  27  -  February  3,  1976)  is  performed  using  a  perfectly 

plastic  model.  Each  calculation  uses  a  different  yield  surface  size 

and/or  shape.  A  teardrop,  triangle  or  diamond  shape  (Figure  2)  is 

4-1  6  ■ 

chosen  and  strength  varies  from  10  Nm  to  10  Nm 


In  addition  a 
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free  drift  calculation  is  performed  to  describe  the  limiting  case  of 
zero  ice  strength. 

The  actual  power  input  to  the  ice  from  the  atmosphere  is  presented 

in  Figure  4.  This  curve  represents  the  average  power  P  over  the 

a 

region  of  interest.  To  calculate  this  variable,  a  set  of  triangles 

connecting  data  buoy  and  manned  camp  locations  is  constructed  as  shown 

T 

in  Figure  1.  The  local  power  p  =  v  T  ,  where  v  is  observed 

a  ~o  ~a  ~o 

station  velocity,  is  then  integrated  over  each  triangle  assuming  linear 

variations  of  p^  within  the  triangle.  Finally,  contributions  from  all 

triangles  are  weighted  by  area  and  summed  to  determine  the  area  average. 

The  power  input  to  the  ice  is  observed  to  be  small  until  the  latter 

2 

half  of  January  29,  rising  to  25  mW/m  late  on  January  30.  The  power 
input  then  drops  to  a  negligible  value  by  the  end  of  January  31  and 
remains  low  through  February  3.  As  the  storm  arrives  (1200  GMT  on 
January  29  through  2100  GMT  on  January  30) ,  the  increase  in  energy  input 
is  caused  largely  by  the  increased  area  of  ice  that  moves.  Winds  and 
ice  speeds  of  the  moving  region  are  less  variable.  But  as  the  storm 
passes,  the  drop  in  power  input  is  caused  by  falling  winds  and  ice 
speeds  over  the  whole  region. 

The  free  drift  power  available  P  ^  is  also  shown  on  Figure  4. 
While  similar  in  shape,  it  is  seen  that  P^^  Is  approximately  twice  as 
large  as  the  observed  value.  The  difference  is  attributed  to  the  effect 
of  ice  stress  since  other  model  results  have  shown  the  boundary  layer 
models  to  be  accurate  at  this  time  (Pritchard,  et  al.,  1977;  Pritchard, 
1977).  The  free  drift  power  available  is  nearly  proportional  to  wind 
speed  cubed  at  higher  winds.  The  increasing  wind  intensity  is  shown 
clearly  in  this  display.  In  fact,  it  is  felt  that  free  drift  power 
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available  is  a  useful  variable  to  determine  which  times  are  of  interest 
for  modeling  ice  behavior.  In  this  display,  storms  are  amplified  to 
show  their  strength  in  a  significant  measure  -  the  power  available  in 
the  region. 

The  three  calculations  titled,  Runs  3B,  C  and  D,  were  reported  by 

Pritchard,  et  al.  (1977).  The  air  stress  used  in  Run  3  is  in  error  on 

January  30  (Pritchard,  et  al. ,  1977).  The  values  input  the  model  are 

too  large  by  about  30  percent,  therefore,  when  model  results  are  compared 

with  observed  values,  they  should  be  scaled  down  accordingly.  The  power 

P^  and  P^  should  be  reduced  to  about  60  percent  of  calculated  values 

for  Run  3  results  on  January  30  because  both  air  stress  and  velocity 

must  be  reduced,  whereas  observed  power  input  should  be  lowered  to 

about  70  percent  because  only  affects  it.  The  velocity  fields  in 

Run  3B  are  judged  to  be  nearly  the  same  as  free-drift.  It  is  seen  that 

the  sequence  of  daily  values  shown  in  Figure  4  for  Run  3B  agrees  to 
2 

within  2  mW/m  with  free  drift  power  available.  This  confirms  the 

4  -1 

conclusion  that  for  yield  strengths  of  10  Nm  the  velocity  field  is 
essentially  free  drift.  On  the  other  hand,  the  velocity  fields  calcu¬ 
lated  in  Run  3C,  using  a  yield  strength  of  10"*  Nm  ^  ,  represent  observed 
velocities  reasonably  well.  The  energy  input  for  Run  3C  (and  the  con¬ 
tinuous  run  using  the  same  model)  exceeds  observed  values  by  about  50 
percent  of  the  peak  observed  value  on  January  30  and  31.  When  scaled 
down  to  correct  for  the  air  stress  error,  the  power  input  is  nearly 
correct . 

As  pointed  out  by  Pritchard  (1977),  the  teardrop  shaped  yield 
surface  incorrectly  simulates  the  uniaxial  opening  of  a  large-scale  lead 
system  because  large  tractions  exist  across  the  open  leads.  It  was 
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shown  that  either  a  triangle  or  diamond  yield  surface  could  correct  this 

deficiency.  When  a  triangle  yield  surface  with  a  strength  p*  *  10^  Nm  * 

(Run  3E)  is  used  to  simulate  ice  behavior  on  January  30,  the  deformation 

and  stress  fields  showed  uniaxial  opening  of  North-South  running  leads 

with  no  traction  across  them.  Furthermore,  a  flaw  lead  is  generated 

about  200  km  offshore,  only  slightly  seaward  of  the  correct  location, 

and  the  velocity  field  is  satisfactory.  The  power  input  from  this  one- 

2 

day  simulation  is  shown  on  Figure  4.  The  power  input  of  12  mW/m  is 
about  30  percent  lower  than  the  observed  value  for  that  day.  This  value 
is  closer  than  previous  calculations,  but,  when  considered  along  with 
the  flaw  lead  position,  indicates  that  yield  strength  is  too  high. 

The  calculation  labeled  Run  3J  gives  results  nearly  identical  to 
Run  3E.  This  occurs  because  the  entire  region  is  either  opening  or 
shearing  with  no  substantial  convergence.  For  these  stretching  states  a 
diamond  with  strength  p*  =  2  x  10^  Nm  *  and  a  triangle  with 
p*  =  1  x  10^  Nm  1  admit  the  same  stress  states  (See  Figure  2).  The 
triangle  is  the  largest  concave  yield  surface  that  is  contained  by  the 
tensile  cutoff  lines  and  a  given  yield  strength  p*  .  Since  an  extreme 
shape  is  unlikely  in  practice,  the  study  uses  the  intermediate  diamond 
shape. 

Three  more  calculations  of  January  30  with  the  diamond  yield  surface 
and  strengths  of  1  x  105  (Run  3H) ,  1.5  x  105  (Run  31)  and  10  x  105  Nm-1 
(Run  3K)  show  how  the  power  input  varies  with  this  parameter  (Figure  4). 
It  is  seen  that  p*  =  1.5  x  10^  Nm  ^  (Run  31)  reproduces  the  observed 
powet  input.  The  results  of  Runs  3H  and  K  provide  a  wider  range  varia¬ 
tion.  From  the  power  input  a  diamond  yield  surface  with 
p*  ■  1,5  x  10~*  Nm  ^  is  best.  The  velocity  and  deformation  fields  for 
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Runs  3H,  I,  J  and  K  show  that  the  flaw  lead  moves  seaward  with  increasing 
strengths  and  that  both  fields  are  accurate  when  strength  is  about  1- 
1.5  x  10^  Nm  1  .  When  the  air  stress  error  is  considered,  the  lower 
strength,  p*  *  1  x  10"*  Nm  ^  is  felt  to  be  best.  The  velocity  field  at 
the  larger  strength  (Run  3J)  is  actually  more  accurate  than  Run  3H 
because  the  eastern  part  of  the  region  is  more  accurate.  When  strength 
is  increased  to  10  x  10“*  Nm  *  (Run  3K)  results  are  similar  to  Run  3D  in 
which  the  flaw  lead  is  far  offshore  and  the  motion  is  dominated  by  buoy 
motions  at  the  boundaries. 

Figure  5  shows  the  power  input  to  the  ocean  P  .  Here  again,  the 

w 

diamond  yield  surface  fits  the  observed  values  well  on  January  30.  The 
best  fit  is  with  p*  =  1  x  10~*  Nm  ^  . 

Tuning  the  Ice  Model 

The  previous  parameter  study  shows  that  the  average  power  input  to 
the  ice  from  the  atmosphere  and  the  power  input  to  the  ocean  P^ 

provide  quantitative  measures  for  verifying  model  performance.  These 
variables  help  describe  how  well  the  velocity  and  deformation  fields  fit 
observed  values.  The  two  calculations  that  most  closely  approximate  the 
power  inputs  are  Runs  3H  and  I. 

Although  the  ideal  plastic  ice  model  just  described  is  used  to 
simulate  a  few  days  of  ice  behavior  accurately,  there  are  cases  where  it 
must  be  inaccurate.  When  a  region  undergoes  large  deformations  (either 
opening  or  closing) ,  changes  in  the  thickness  distribution  (production 
or  depletion  of  open  water  and  thin  ice)  must  change  the  yield  strength 
locally.  If  this  softening  or  hardening  is  ignored,  as  in  an  ideal 
plastic  model,  then  subsequent  response  can  be  in  serious  error.  With 


this  in  mind,  the  remainder  of  this  work  uses  results  of  a  model  that 


allows  hardening/ softening  behavior.  Strength  changes  depend  on  thick¬ 
ness  distribution  changes.  This  concept  is  an  important  physical 
feature  of  the  AIDJEX  model  (Coon  et  al.,  1974).  Recently,  a  modified 
formulation  has  allowed  this  hardening  behavior  to  be  included  in  a 
model  that  has  strength  on  the  order  of  10^  Nm  ^  (Coon,  et  al . ,  1978). 
Such  strengths  were  not  compatible  with  physical  definitions  of  material 
constants  in  the  original  AIDJEX  model.  This  final  version  of  the 
AIDJEX  model  includes  a  shear  energy  sink  that  treats  energy  dissipation 
without  redistribution  of  ice.  This  generalization  and  a  redefinition 
of  material  constants  increase  the  initial  yield  strength  to  1  x  10"*  Nm 
and  also  prohibit  softening  in  many  shearing  situations  where  the 
previous  model  allowed  it.  In  addition,  thickness  distribution  has  been 
modified  to  represent  the  latest  keel  depth  distributions  measured  by 
the  U.S.S.  Gurnard  upward-looking  sonar  in  April  1976  (Wadhams,  1978). 
These  values  are  extrapolated  into  the  January-February  1976  time  frame 
using  the  AIDJEX  redistribution  model.  Our  estimate  of  actual  ice 
conditions  were  presented  by  Coon  et  al.  (1978). 
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Understanding  the  Winter  Energy  Budget  in  the  Beaufort  Sea 

In  this  section  the  energy  budget  for  the  hardening,  plastic  ice 

model  with  tuned  parameters  is  studied.  Particular  interest  is  given 

to  modification  of  the  rate  of  energy  transmitted  into  the  ocean  and 

to  the  horizontal  transfer  of  energy  by  the  ice  stress. 

As  discussed  in  the  Section  on  test  conditions,  two  storms  passed 

through  the  Beaufort  Sea  during  the  17  day  interval  (Figure  6) .  The 

first  storm  reached  a  maximum  during  January  30,  and  the  second  storm 

reached  a  maximum  during  February  7  and  8.  Time  histories  of  P  and 

P^  provide  a  vivid  picture  of  the  intensity  of  these  storms. 

The  power  input  from  the  atmosphere  to  the  ice  P^  is  determined 

from  the  winds  and  motions  for  the  triangles  in  Figure  1.  The  power  to 

the  ocean  P  is  determined  from  the  motion  and  ocean  stress  model, 
w 

Figure  7  shows  the  P^  and  P^  calculated  from  a  run  (Run  5C)  of  the 

tuned  hardening/softening  model.  Comparing  Figure  6  and  7  it  can  be 

seen  that  these  two  terms  are  approximated  very  well.  Figure  7  also 

shows  a  curve  for  P^  +  P^  when  P^  is  the  rate  at  which  energy  comes 

in  to  or  out  of  the  region  from  the  ice  to  the  North  and  West  of  the 

region.  It  is  the  combination  of  P  +  P  that  is  available  to  deform 

a  i 

the  ice  and  go  into  the  ocean.  On  January  30,  P  +  Pf  is  less  than 

a  x 

P  and  energy  is  leaving  the  region.  During  the  storm  of  February  7 
a 

and  8  energy  is  coming  into  the  region.  On  February  6  P  +  Pf  is 

a  I 

greater  than  P  ^  and  therefore  more  energy  is  coming  into  the  region 

than  would  be  available  in  the  region  for  a  free  drift  (no  strength)  ice 

condition.  The  difference  between  curves  P  +  P-  and  P  is  the 

a  f  w 

power  that  is  dissipated  by  the  ice.  We  select  January  30  and  February  7 
for  further  study  of  the  spatial  distribution  of  the  various  power 
density  measures  in  this  work. 
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The  power  available  from  the  atmosphere  is  expressed  in  terms  of 
the  free  drift  power  density,  and  may  be  seen  in  Figure  8  for  January  30, 

It  is  seen  that  the  maximum  power  available  is  located  approximately 
300  km  north  of  the  barrier  islands  off  the  north  slope  of  Alaska,  and 
small  amounts  of  energy  are  available  locally  along  the  north  slope. 

The  circular  contours  represent  constant  levels  of  wind  speed,  and  are 
approximately  proportional  to  the  cube  of  the  wind  speed.  The  wind  field  is 
not  directly  influenced  by  presence  of  the  land  mass.  Air  stress  vectors 

are  also  shown  on  Figure  8.  The  maximum  power  density  available  is 

2  2 
60  mW/m  ,  the  power  density  near  the  shore  ranges  from  10-20  mW/m 

The  actual  power  input  from  the  atmosphere  is  shown  in  Figure  9 

where  wind  and  buoy  motions  are  used  in  the  calculation.  Data  points 

at  each  buoy  are  shown.  There  is  no  power  input  along  the  north  slope 

2 

because  the  ice  is  stationary.  The  largest  input  power  is  55  mW/m 

The  location  of  the  fast  ice  region  is  outlined  by  the  dashed  line. 

Small  variations  shoreward  of  the  dashed  line  are  ignored  in  drawing 

contours.  When  the  same  air  stress  field  is  used  as  input  to  the  model 

simulation,  the  resulting  velocity  field  provides  the  actual  power  input 

to  the  ice  as  shown  in  Figure  10  for  January  30.  It  is  seen  that  the 

same  elliptical  pattern  of  contours  remains  and  the  observed  field  is 

approximated  reasonably  well.  The  maximum  power  density  input  to  the 

2 

ice  at  this  time  is  30  tnW/m  with  a  minimum  along  shore  similar  to  that 

2 

available  in  the  atmosphere  and  on  the  order  of  10  mW/m 

The  measured  power  transmitted  to  the  ocean  on  January  30  Is  shown 
in  Figure  11 .  This  power  is  calculated  using  buoy  motion  and  the  water 
stress  model.  The  amount  of  power  p^  calculated  by  the  ice  model  that 
is  transferred  into  the  ocean  from  the  ice  cover  is  shown  in  Figure  12 


2 

for  January  30.  This  compares  within  3  mW/m  with  the  measured  values 

of  Figure  11.  The  magnitudes  are  dramatically  smaller  than  that  input 

into  the  ice  from  the  atmosphere.  Maximum  values  input  are  in  the  order 
2 

of  10  mW/m  with  almost  the  entire  study  area  receiving  less  than 
2 

5  mW/m  .  The  power  input  to  the  ocean  from  the  ice  is  proportional  to 
the  cube  of  the  ice  speed  if  we  ignore  geostrophic  ocean  currents.  The 
superposed  velocity  field  indicates  larger  speeds  in  the  western  portion 
of  the  basin  and  smaller  speeds  in  the  eastern  region  that  has  just 
begun  to  move  and  these  variations  are  reflected  in  the  contours  of 
Figure  11. 

The  difference  between  power  input  from  the  atmosphere  and  that 

transmitted  into  the  ocean  is  p^  =  pa  *  p^  .  This  power  represents  the 

effect  of  ice  stress.  On  January  30  it  is  shown  in  Figure  13.  Since 

the  term  is  much  larger  than  p^  Figure  13  tends  to  look  very 

similar  to  the  power  from  the  atmosphere  (Figure  10) .  It  is  seen  that 

the  largest  power  input  t<i  the  ice  is  in  the  middle  of  the  region  and  is 

2 

in  the  order  of  30  mW/m  .  Small  values  are  input  locally  near  the 
edges  of  the  basin. 

This  power  input  to  the  ice  p^  may  be  decomposed  into  the  power 
dissipated  by  deformation  p^  and  the  divergence  of  the  stress  flux 
p^  .  The  power  from  the  divergence  of  the  stress  flux  is  shown  in 
Figure  14.  In  this  figure  there  are  both  negative  and  positive  contours. 
A  negative  value  shows  that  energy  is  leaving  the  point  and  a  positive 
value  shows  that  energy  is  coming  into  the  point  through  the  ice  cover. 
The  area  integral  of  the  power  from  the  divergence  of  the  stress  flux  on 
January  30  must  be  negative.  This  is  shown  in  Figure  7,  which  presents 
the  average  power  over  the  region.  At  this  time  energy  is  leaving  the 


region  at  the  west  boundary  with  the  stress  in  compression  and  the 
motion  of  the  ice  to  the  west.  The  high  positive  values  of  contours 
along  the  north  slope  show  that  a  large  amount  of  energy  is  coming  into 
that  area. 

The  stress  power  or  power  dissipated  by  deformation  p.  on  Janu- 

d 

ary  30  is  shown  in  Figure  15.  In  contrast  to  the  power  input  locally  to 

the  ice  by  the  difference  between  atmosphere  and  ocean  terms ,  the  power 

dissipated  by  the  ice  cover  is  concentrated  along  the  shorelines.  In 

the  central  Arctic  regions,  it  is  seen  that  dissipation  by  the  ice  cover 

is  zero  (this  term  cannot  be  less  than  zero)  whereas  along  the  shore- 

2 

lines,  a  maximum  on  the  order  of  100  mW/m  is  obtained.  These  large 
values  occur  when  large  deformations  along  shorelines  and  large  stresses 
exist  simultaneously.  In  Figure  15  dissipation  in  the  vicinity  of  75° 
North  and  130°  West  is  introduced  by  spurious  motion  of  the  boundary 
buoy  and  should  be  ignored.  However,  the  concentration  of  large  amounts 
of  dissipation  along  the  north  slope  of  Alaska  is  physically  real.  The 
band  of  shearing  that  causes  this  large  energy  dissipation  is  seen  to  be 
approximately  150  km  wide  and  covers  approximately  three  cells  in  the 
numerical  grid.  For  the  calculation  discussed  in  this  work  the  maximum 
dissipation  occurs  at  the  shoreline.  Qualitatively  this  result  is  quite 
acceptable.  However,  other  simulations  (Pritchard,  1977)  have  shown 
that  the  flaw  lead  existed  at  this  time  approximately  150  km  offshore. 
Quantitatively,  therefore,  one  should  expect  that  there  is  very  little 
energy  dissipated  shoreward  of  this  150  km  line,  and  the  concentrations 
seen  in  Figure  15  should  be  moved  out  to  the  location  of  the  flaw  lead. 
It  is  known  from  previous  parameter  studies  that  this  offshore  flaw  lead 
is  modeled  accurately  when  strengths  are  increased  to  the  order  of 


1.5  x  103  Nm  .  Those  higher  strength  values  are  also  consistent  with 

results  of  the  parameter  study  discussed  earlier  in  this  work.  It  is 

seen  that  the  power  available  from  the  atmosphere  p  and  input  to  the 

a 

ice  p^  over  the  entire  Beaufort  Sea  (Figures  10  and  13)  is  transferred 
by  the  stress  in  the  ice  cover  and  dissipated  over  the  Continental 
Shelf.  This  physical  process  of  a  horizontal  energy  transfer  is  dif¬ 
ferent  from  open  oceans  and  appears  to  be  an  important  process  in  ice 
covered  oceans.  It  is  a  mechanism  that  controls  the  dissipation  of 
marine  energy  over  the  shelf  regions  in  polar  oceans.  It  should  be 
pointed  out  that  along  the  eastern  shore  of  the  region,  near  Banks 
Island,  the  power  dissipated  by  the  ice  cover  is  negligible  even  though 
deformations  in  that  region  are  large.  The  opening  of  a  lead  along  the 
land  is  occurring.  The  small  power  dissipation  occurs  because  stresses 
are  negligible  in  that  opening  region. 

If  the  generation  of  acoustic  energy  by  ice  ridging  is  of  concern, 
the  power  dissipation  contours  can  be  used  to  determine  where  energy  is 
being  dissipated  by  the  ice  cover.  It  is  easy  to  imagine  that  the 
energy  dissipated  by  deformation  will  create  noise  and  much  of  this 
noise  will  be  transmitted  into  the  ocean  and  propagated  from  that  point. 
Operators  in  the  Beaufort  Sea  should  well  suspect  that  on  January  30 
conditions  over  the  Continental  Shelf  along  north  slope  of  Alaska  would 
be  quite  noisy.  By  the  same  token,  one  should  expect  relatively  little 
contribution  to  the  noise  level  from  the  area  near  Banks  Island. 

The  conditions  described  in  Figures  8-15  are  easterly  winds  and  ice 
velocities  toward  the  west.  Large  shearing  deformations  along  the  north 
slope  of  Alaska  and  uniaxial  opening  along  the  shoreline  of  Banks  Island 
typically  occur  under  these  conditions.  In  the  next  sequence  of  figures. 


conditions  where  ice  drift  and  winds  are  towards  the  McKenzie  delta  and 

Banks  Island  are  presented.  These  create  convergence  in  the  Canadian 

Beaufort  Sea.  At  the  same  time  shearing  occurs  along  the  Alaskan  north 

slope  but  the  shearing  is  in  the  opposite  sense  compared  with 

previous  days.  The  power  available  from  the  atmosphere  at  this  time  is 

shown  in  Figure  16  for  February  7.  For  this  case,  contour  levels  run 

approximately  from  the  southwest  to  the  northeast,  with  maximum  winds  to 

the  north  of  Barter  Island.  The  power  available  at  this  time  is  a 

2 

maximum  of  approximately  30  mW/m  .  The  average  power  available  over 

2 

the  entire  region  of  interest  is  approximately  27  mW/m 

The  observed  p  is  shown  in  Figure  17  and  the  p  from  the  model 
a  a 

calculation  are  shown  in  Figure  18.  The  power  input  to  the  ice  from  the 

atmosphere  on  February  7  is  affected  significantly  by  the  presence  of 

the  land  mass.  (Compare  Figures  16  and  17.)  Contours  are  seen  to 

generally  follow  the  shape  of  the  north  slope  of  Alaska  and  Canada. 

Maximum  values  input  to  the  ice  from  the  atmosphere  according  to  this 

2 

simulation  are  on  the  order  of  20  mW/m  .  Since  the  ice  is  assumed  to 

be  at  rest  at  the  shore,  one  must  always  expect  zero  input  at  this  point 

because  of  the  zero  velocity.  The  exception  to  this  case  arises,  of 

course,  when  a  velocity  discontinuity  appears  at  the  shoreline,  but  the 

numerical  solution  scheme  does  not  represent  such  discontinuities. 

The  measured  p  is  shown  in  Figure  19  and  the  p  from  the  model 

calculation  is  shown  in  Figure  20.  The  power  input  to  the  ocean  on 

February  7  is  small  throughout  the  region  of  interest.  Except  at  the 

2 

extreme  northwest  boundary,  power  throughout  the  area  is  5  mW/m  or 
less.  The  small  values  occur  because  the  velocity  throughout  the  region 
is  small  as  seen  in  the  superposed  Ice  drift  vectors  (Figure  20) .  The 
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velocity  field  is  small  in  magnitude  at  this  time  because  onshore  winds 

during  the  last  day  have  caused  convergence  of  the  ice  cover  along  the 

eastern  shore  and  the  material  has  hardened.  Hardening  of  the  ice  model 

allows  an  increase  in  stress  levels  which  prevents  further  motion. 

The  effect  of  the  ice  for  February  7  is  shown  in  Figure  21  with 
2 

contours  from  5-15  mW/m  .  The  effect  of  the  ice  is  uniform  over  the 
region.  However,  when  is  decomposed  into  and  p^  this  view 

changes.  The  power  from  the  divergence  of  the  stress  flux  on  February  7 
is  shown  in  Figure  22.  Here  again,  both  negative  and  positive  contours 
occur.  However,  there  is  a  net  energy  input  to  the  region  as  can  be 
seen  from  Figure  7.  The  areas  along  the  north  slope  and  along  the 
Canadian  Archipelago  show  that  large  amounts  of  energy  are  being  trans¬ 
ferred  to  these  areas. 

Power  dissipated  by  the  ice  cover  on  February  7  is  presented  in 
Figure  23.  As  usual,  it  is  concentrated  around  the  shorelines.  Along 
the  north  slope  of  Alaska  a  strong  shearing  occurs  as  the  ice  moves 

2 

toward  the  east  with  local  power  dissipation  on  the  order  of  80  mW/m 

These  values  fall  rapidly  within  the  first  75  km  from  shore.  Along 

Banks  Island  and  Tuktoyaktuk  Peninsula  the  convergence  again  creates 

large  values  of  dissipation.  Maximum  local  values  occur  across  Amundsen 

2 

Gulf  and  reach  a  magnitude  of  50  mW/m  .  In  this  region  of  uniaxial 
convergence,  contours  are  spread  wider  spatially  than  they  are  in  the 
region  of  shear  along  the  Alaskan  north  slope.  This  spreading  is  caused 
by  hardening  that  occurs  in  compression.  As  the  material  hardens, 
dissipation  occurs  throughout  the  region  of  hardening.  In  contrast, 
shearing  occurs  with  little  change  in  yield  strength,  and  therefore 
deformations  are  concentrated  either  along  discontinuities  or  narrow 
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Conclusion 

The  mechanical  energy  variables  considered  in  this  work  provide 
useful  measures  of  comparison  between  modeled  and  observed  sea  ice 
behavior.  The  power  input  to  the  ice  cover  from  the  atmosphere  is 
useful  for  determining  when  strong  storms  or  interesting  extreme  events 
occur. 

If  barometric  pressure  fields  (or  wind  fields)  and  ice  motion  data 
are  available,  the  free  drift  power  available  from  the  atmosphere,  the 
actual  power  input  to  the  ice  from  the  atmosphere  and  the  power  input  to 
the  ocean  can  be  found.  Although  useful,  these  variables  are  inadequate 
for  estimating  other  important  terms  in  the  mechanical  energy  budget. 

To  determine  the  mechanical  energy  budget  completely  requires  an 
ice  model  or  constitutive  law.  With  this  model,  knowledge  of  the  ice 
stress  allows  stress  flux  divergence  and  energy  dissipation  by  the  ice 
to  be  evaluated.  The  stress  flux  divergence  is  needed  to  determine  how 
energy  is  transferred  horizontally  through  the  ice  cover  and  explains 
the  fact  that  energy  is  dissipated  primarily  around  the  margins  of  the 
ice-covered  oceans.  These  processes  should  be  helpful  in  interpreting 
how  coastal  erosion  is  altered  by  the  ice  cover  and  how  acoustic  energy 
is  generated  in  these  areas, 

A  comparison  of  energy  input  to  the  ice  from  the  atmosphere  and  of 
the  energy  transferred  to  the  ocean  has  allowed  a  plastic  sea  ice  model 
to  be  tuned.  We  have  first  of  all  shown  that  these  energy  measures  are 
useful  for  model  tuning  and  have  then  used  the  results.  A  diamond  shape 
yield  surface  is  seen  to  be  better  than  the  teardrop  used  previously.  A 
yield  strength  of  1  x  10^  Nm  ^  is  judged  best.  However,  this  is  slightly 
different  from  the  value  of  1.5  x  10^  Nm  ^  that  would  be  chosen  by 


> 


f 


considering  only  the  velocity  field.  Such  differences  point  out  the 
importance  of  comparing  all  quantities  of  interest  rather  than  just 
one.  It  is  now  felt  that  yield  strength  and  yield  surface  shape  are 
defined  well  enough  that  further  work  must  consider  many  more  physical 
processes  simultaneously.  These  include  anisotropy,  tensile  strength, 
ice  age,  floe  size  and  scale  effects. 

Results  of  the  model  calculation  during  the  seventeen  day  test 
period  show  that  a  variety  of  conditions  can  arise.  We  have  found  times 
when  energy  is  increased  in  the  Beaufort  Sea  from  distant  regions  and 
other  times  when  it  is  decreased.  At  times  these  boundary  effects  are 
large  enough  so  that  power  input  even  exceeds  maximum  power  available 
from  the  atmosphere  over  the  Beaufort  Sea. 
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Figure  1. 


Beaufort  Sea  Region  Subdivided  into  Triangles  by  Manned 


Camps  and  Data  Buoys. 

Figure  2.  Schematic  of  Important  Terms  Contributing  to  Mechanical 

Energy  Budget. 

Figure  3,  Yield  Surface  Shapes  Considered  in  Parameter  Study. 

Figure  A.  Comparison  of  Modeled  Average  Atmospheric  Power  to  Ice  in 

Beaufort  Sea  for  a  Range  of  Strength  and  Yield  Surfaces. 
Figure  3.  Comparison  of  Modeled  Power  to  Ocean  in  Beaufort  Sea  for 

a  Range  of  Strength  and  Yield  Surfaces. 

Figure  6.  Observed  Averaged  Atmosphere  Power  to  Ice  and  to  Ocean  in 

Beaufort  Sea  During  Seventeen  Day  Test  Period. 

Figure  7.  Average  Power  Variables  Calculated  from  the  Best  Model 

During  Seventeen  Day  Test  Period. 

Figure  8.  Free  Drift  Power  Available  p  ^  on  January  30.  All 

2 

contour  levels  in  mW/m  .  Scale  vector  is  for  air 
stress  vectors. 

Figure  9.  Observed  Atmospheric  Power  to  Ice  p  on  January  30. 

2 

All  contour  levels  in  mW/m 

Figure  10.  Modeled  Atmospheric  Power  to  Ice  p  on  January  30.  All 

a 

2 

contour  levels  in  mW/m 

Figure  11.  Observed  Power  to  Ocean  p  on  January  30.  All  contour 

2 

levels  in  mW/m 

Figure  12.  Modeled  Power  to  Ocean  p^  on  January  30.  All  contour 

2 

levels  in  mW/m  .  Scale  vector  is  for  ice  velocity 


vectors. 
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Figure  13. 


Figure  14. 


Figure  15. 


Figure  16. 


Figure  17. 


Figure  18. 


Figure  19. 


Figure  20. 


Figure  21. 


Figure  22. 
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Free  Drift  Power  Available  Pafd  on  February  7.  All 

2 

contour  levels  in  mW/m  .  Scale  vector  is  for  air  stress 
vectors. 

Observed  Atmospheric  Power  to  Ice  p^  on  February  7. 

2 

All  contour  levels  in  mW/m 

Modeled  Atmospheric  Power  to  Ice  p  on  February  7.  All 

2 

contour  levels  in  mW/m 

Observed  Power  to  Ocean  p^  on  February  7.  All  contour 
2 

levels  in  mW/m 

Modeled  Power  to  Ocean  p^  on  February  7.  All  contour 
2 

levels  in  mW/m  .  Scale  vector  is  for  ice  velocity 
vectors. 

Effect  of  Ice  Stress  p^  on  February  7.  All  contour 
2 

levels  in  mW/m 

Stress  Flux  Divergence  p^  on  February  7.  All  contour 
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Figure  23. 


Beaufort  Sea  Region  Subdivided  into  Triangles  by  Manned  Camps  and  Data  Buoys. 
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Comparison  of  Modeled  Average  Atmospheric  Power  to  Ice  in  Beaufort  Sea 
for  a  Range  of  Strength  and  Yield  Surfaces. 
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Figure  5.  Comparison  of  Modeled  Power  to  Ocean  in  Beaufort  Sea  for  a  Range  of 
Strength  and  Yield  Surfaces. 


•gure  6.  Observed  Averaged  Atmosphere  Power  to  Ice  and  to  Ocean  in  Beaufort  Sea 
During  Seventeen  Day  Test  Period. 


Figure  7.  Average  Power  Variables  Calculated  from  the  Best  Model  During  Seventeen 
Day  Test  Period. 


Figure  8.  Free  Drift  Power  Available  pafd  on  January  30.  All  Contour  Levels  in 
mW/m^ .  Scale  Vector  is  for  Air  Stress  Vectors. 


Figure  9.  Observed  Atmospheric  Power  to  Ice  pa  on  January  30.  All  Contour  Levels 
in  mW/m^ 


Figure  10.  Modeled  Atmospheric  Power  to  Ice  pa  on  January  30.  All  Contour  Levels 
in  mW/m^. 


Figure  12.  Modeled  Power  to  Ocean  pw  on  January  30.  All  Contour  Levels  in  mW/nv 
Scale  Vector  is  for  Ice  Velocity  Vectors. 
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Figuare  14.  Stress  Flux  Divergence  pf  on  January  30.  All  Contour  Levels  in  mW/m^ . 


Flguare  15.  Power  Dissipated  by  Ice  on  January  30.  All  Contour  Levels  in  mW/m^ . 
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Modeled  Atmospheric  Power  to  Ice  pa  on  February  7.  All  Contour  Levels  in  mW/m  . 


Figure  19.  Observed  Power  to  Ocean  pw  on  February  7.  All  Contour  Levels  in  mW/m^ 


Scale  Vector  is  for  Ice  Velocity  Vectors. 


Figure  21.  Effect  of  Ice  Stress  p*,  on  February  7.  All  Contour  Levels  in  mW/m^ 


Figure  22.  Stress  Flux  Divergence  pj  on  February  7.  All  Contour  Levels  in  mW/m^ 


Figure  23.  Power  Dissipated  by  Ice  p,j  on  February  7.  All  Contour  Levels  in  mW/m^ 
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